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Abstract 

We describe a midi-superspace quantization scheme for generic single 
horizon black holes in which only the spatial diffeomorphisms are fixed. The 
remaining Hamiltonian constraint yields an infinite set of decoupled eigen- 
value equations: one at each spatial point. The corresponding operator at 
each point is the product of the outgoing and ingoing null convergences, and 
describes the scale invariant quantum mechanics of a particle moving in an 
attractive potential. The variable X that is analoguous to particle po- 

sition is the square root of the conformal mode of the metric. We quantize 
the theory via Bohr quantization, which by construction turns the Hamil- 
tonian constraint eigenvalue equation into a finite difference equation. The 
resulting spectrum gives rise to a discrete spatial topology exterior to the 
horizon. The spectrum approaches the continuum in the asymptotic region. 
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1 Introduction 



It is commonly believed that the puzzling thermodynamic properties of black 
holes have their origin in the quantum behaviour of the gravitational field. 
The two leading candidates for a consistent quantum theory of gravity each 
have their own picture of the structure that underlies spacetime geometry. In 
the case of string theory, the fundamental microscopic objects are supersym- 
metric strings and branes, whereas in loop quantum gravity the fundamental 
structures are spin networks that generate a quantum phase space whose co- 
ordinates are holonomies on the spin networks and densitized spatial triads. 
Although both theories have made progress within their respective frame- 
works^ , it is safe to say that a great deal remains to be learned, particularly 
about the quantum gravitational origins of the thermodynamic properties of 
black holes. 

One striking aspect of black hole thermodynamics is the apparent uni- 
versality: the key features seem not to depend on the specific action from 
which they are derived. All that is required is the presence of an event hori- 
zon in a diffeomorphism invariant theory. It therefore makes sense to try to 
understand the quantum behaviour of black holes, and the resulting ther- 
modynamics not from the point of view of a specific microscopic underlying 
theory, but instead in terms of the things that they have in common (at the 
semi-classical level at least), namely their geometrodynamic properties. 

Recently Husain and Winkler have initiated a program |2j designed to 
address directly the quantum geometrodynamics of black holes. Their view- 
point is that returning to geometrodynamical variables opens up new pos- 
sibilies for resolving the Hamiltonian constraint. Of course, this cannot be 
accomplished by a return to the original formulation, wherein the quantum 
phase space was parametrized by the spatial metric and its conjugate momen- 
tum. Instead, adapting a quantization method used recently in the context 
of loop quantum gravity(see 0), the latter is replaced by its exponentiated 
version. This necessitates the construction of a quantum phase space with a 
discrete topology, as is the case with 'conventional' loop quantum gravity. 

At the present time, this program has only been fully implemented for 
'reduced' systems, that is, for the case of dynamics of homogeneous and 
isotropic universes, or where there is spherical symmetry. Since a good part 

"'^See for example in which it is shown that loop quantum cosmology is capable of 
resolving the big bang singularity 
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of gravitational physics deals with these two cases, there have already been 
interesting results regarding singularity resolution and in understanding the 
quantum definition of horizons |2]. The first steps have also been taken 
towards implementation of the program in a partially reduced theory with 
matter 

In this paper, we will, first of all, broaden the applicability of this program 
by applying it to generic single horizon black holes using the formalism of two 
dimensional dilaton gravity models. This extends the program to spherically 
symmetric gravity in arbitrary dimensions, as well as to models of interest 
in string theory, such as the CGHS model of two dimensional gravity 

Second of all, we will use a method of gauge fixing distinct from |1] that 
turns out to have interesting properties as well as potential implications be- 
yond the Bohr quantization program and the vacuum case considered here. 
In particular, we look at a midi-superspace model obtained by partially fix- 
ing the spatial diffeomorphisms, so that in effect, the spatial coordinate is a 
specified function of the dilaton field. Our partial gauge fixing is very nat- 
ural from the geometrodynamics viewpoint and allows one to examine the 
evolution of spatial slices that extend from the singularity across the horizon 
and off to infinity. Remarkably, the resulting Hamiltonian operator reduces 
to the product of the ingoing and outgoing null convergences at each spatial 
point. Moreover, the Hamiltonian at each point corresponds to that of a 
particle moving in a potential. The potential is scale invariant 

and its quantization has been studied extensively in the past [HI El , both for 
its interesting mathematical properties and for its physical relevance. Most 
significantly, it arises in the near horizon dynamics of black holes where the 
conformal symmetry has been used to account for the entropy associated with 
the horizonPlllUj. (See also jH] for an analysis of the connection between this 
potential and black hole thermodynamics.) 

We recognize the relative triviality of vacuum dilaton gravity- there are 
no propogating degrees of freedom. Nevertheless, we will show that the im- 
plementation of the Bohr quantization program in this context allows us to 
solve the Hamiltonian constraint, and sets the stage for our consideration 
of quantum evolution of systems containing gravity and propagating mat- 
ter. Moreover, if indeed diffeomorphisms near the boundary do play a role in 
black hole entropy, as suggested by Carlip^ and others, it is of great interest 
to able to quantize these modes as well as the physical ones. Finally, and per- 
haps most significantly, our quantization scheme has interesting consequences 
even for black hole vacuum spacetime: it induces a discrete spatial topology 
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exterior to the horizon that approaches the continuum at large distances. 



2 Classical Theory 



We start from the action 




) 



(1) 



In the above G is the two dimensional Newton's constant, which has units 
of 1/h, while / is an arbitrary length scale that we take to be the Planck 
scale for the theory. This model has been extensively studied in the past 
using a variety of formalisms (see jTTj for a recent review). This theory is 
exactly solvable both classically and quantum mechanically. For details of 
the analysis in terms of geometrodynamic variables, see ^2] ■ 

This action describes, for different choices of dilaton potential, a very large 
class of theories. Spherically symmetric Einstein gravity in four spacetime 
dimensions corresponds to a dilaton potential V oc 1/ \/0 and on shell = 
is the area of symmetric spheres of radius r. For the more general case of d- 
dimensional Einstein gravity with — 2 dimensional spherical symmetry, V oc 
0-i/(d-2)^ C)C is the area of a unit n-sphere and G = SirG'^'^'^ / , 
where fi'-'^-' is the volume of the unit n-sphere. One can without loss of 
generality take / = (/iG'-'^))^/'-'^"^) = Ipi to be the higher dimensional Planck 
length. 

The vacuum theory admits a one parameter family of solutions with at 
least one Killing vector. This solution can be represented in Schwarzschild 
like coordinates x = 1(f) as: 



where = J''' d(j)V{(f)) and M is the ADM mass. For the present purposes 
it is convenient to assume that is a monotonic function of (p, in which 
case each vacuum solution contains a horizon located at (ph = j~^{2GMl) 

The canonical analysis follows that in [T^. The metric is written in 
modified ADM form: 



We are assuming that the 2D spacetime manifold M2 = i? x S, where S is a 
spacelike slice. As usual the lapse and shift functions a and A^, respectively. 



ds^ = - 2GMl)dt^ + (j(0) - 2GMl)-^dx\ 




(3) 
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are Lagrange multipliers, while the momenta conjugate to the fields and p 
defined above are: 



{Np' + N'-p), 



and the hamiltonian is 



i7 = /(a^+AA^) + surface terms. 



(4) 



(5) 



Thus the Lagange multipliers o", enforce the secondary constraints (respec- 
tively the Hamiltonian and diffeomorphism constraints): 



GG : - -GI[pGY[^ + <p" ■ 



hp 



I 
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The momenta canonically conjugate to a, N are primary constraints. 
An improved Hamiltonian constraint is 

g -le-^P {(f)'g + GH^JT) . 
It turns out that Q — H' is a, total spatial derivative , with 



{g'uI - {cP'f) + 



jM 

P 



(6) 
(7) 



(8) 



(9) 



on the constraint surface. It follows that the Hamiltonian constraint is equiv- 
alent to 

= M = constant. (10) 

where M is the energy, or black hole mass. 

We now perform a canonical transformation: 



so that 



-XP' + 



2G 



X2 



+ 



P 



(11) 

(12) 

(13) 
(14) 
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3 The Classical Null Expansion Observable 



We now verify an important property of the improved Hamiltonian constraint 
operator (JU, namely its relationship to the expansion of null rays. Given 
that we are working in two spacetime dimensions, the definition of expansion 
requires a bit of care. Guidance is provided by the fact that one can consider 
2D dilaton gravity as dimensionally reduced Einstein gravity in D = n + 2 
spacetime dimensions. In this case the dilaton has the interpretation as the 
invariant area of the n— sphere spanned by the generators of spatial rotations, 
i.e. the area of a sphere of fixed radius r. Even in the case where the 2D 
theory is considered fundamental, the dilaton provides the only scalar field 
whose rate of change along null vectors has an invariant meaning. With this 
motivation we define the expansion as the fractional change in the value 
of the dilaton along ingoing and outgoing null vectors, respectively: 

e<*' = ^/w. (15) 

where the vectors 

=n^± s,, (16) 

are ingoing and outgoing null vectors defined in terms of the future pointing 
normal = e^[(T, 0] of each spacelike slice defined by the metric Q, and the 
embedded inward pointing normal = eP[—N, —1]. 



A straightforward calculation then gives: 



e(±) 




(17) 



As in Husain and Winkler ^3|; the expansions provide phase space func- 
tions that allow one to test for the presence of horizons. We now see that the 
Hamiltonian operator can be written in terms of the expansion observables: 



n oc 



e+0- 
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4 Gauge Fixing 



It is well-known that 2D dilaton gravity has no propogating degrees of free- 
dom. This follows immediately from the fact that the Hamiltonian is a linear 
combination of the secondary constraints JF, Q. Thus we in can principle re- 
duce our system to a quantum mechanical system by gauge fixing. Here, we 
fix only the freedom associated with spatial diffeomorphisms by choosing 

:= 0' - = 0. (19) 

which implicitly specifies the spatial coordinate in terms of the invariant dila- 
ton. Recall that is the invariant area of a symmetric n-sphere in spherically 
symmetric gravity in d=n+2 dimensions. At this stage any monotonic func- 
tion f{(j)) would do but we would like to describe the quantum mechanics 
of a black hole spacetime inside, outside and in the neighbourhood of the 
horizon. We therefore make the choice = j{4>), which is consistent with 
the 2D dilaton gravity anologue of Painleve-Gullstrand coordinates, in which 
the generic metric takes the form: 

ds'' = (^-dt' + (dx + ^J^^dtf^ , (20) 

where, from (fTUj) . dx = ld(j)/j{(t>). 

When we smear the diffeomorphism constraint with a test function that 
vanishes on the boundaries: 

X(A) = / dxx{x)X{x), (21) 
and compute its Poisson bracket with the constraint, we find 

{xW,Hy)} = I dxrm [\' + . (22) 

Given the boundary conditions on \{x) this bracket is invertible. Hence, the 
gauge choice above is acceptable, in the Dirac sense. 

The consistency condition that preserves the gauge fixing condition in 
time is: 

aGIip = N(j)', (23) 
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and one can now eliminate from the dynamical equations by solving the 
diffeomorphism constraint. The partially reduced Hamiltonian is: 

H = - J^dxaW + anloundary ■ (24) 

It is easy to show that the Dirac bracket of X with P reduces to the Poisson 
bracket of these variables so that the equations of motion for the remaining 
fields can be derived directly from ()24|) using the ordinary Poisson brackets. 
In this gauge, the theory is determined entirely by the function Ti, which on 
the constraint surface is a constant of motion. The Hamiltonian constraint 
()1U|) now takes the simple form: 

0{X, P) ^ {GlPf - ^ - {2GMI - m) = 0. (25) 

Since we have fixed the spatial diffeomorphisms, the value of at each 
spatial point can be considered a c-number. Thus only P and X are to 
be quantized. In this gauge, there is no coupling between different spatial 
points so that, remarkably, the midi-superspace model has decoupled into 
an infinite set of quantum mechanical models. Moreover, for each value of 
0, the Hamiltonian in is that of a particle moving in a binding 
potential with total "energy" 

e = 2GMl-j{(P), (26) 

determined by the black hole mass and the invariant location along the spatial 
slice as determined by the value of (p. The total "energy" is negative in 
the asymptotic region outside the horizon, so that one is looking for bound 
states in this region, whereas the energy is positive in the interior where one 
is looking for "scattering states". At the classical horizon location j(0) = 
2GMI, the eigenvalue is zero. 

The l/X^ potential arises in a variety of contexts, including the quan- 
tum Hall effect and a field theory formulation of nucleon interactions [T3] . 
More importantly for the present context it arises in the near horizon dynam- 
ics of black holes where the conformal symmetry has been used to account 
for the entropy associated with the horizon [1^ llUj. 

An important property of the inverse square potential Hamiltonian (j^^ 
is its scale invariance. As a result, it possesses three symmetry generators 
that span an algebra isomorphic to the conformal group 5*0(1, 2). The quan- 
tum mechanics of this potential has been the subject of much work, in part 
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because the conformal symmetry is necessarily broken at the quantum level. 
An elegant way of understanding the anomaly was given by EstevejTTj: the 
other generators of the conformal group do not preserve the domain of self 
adjointness of the Hamiltonian, which resulting non-self-adjointness leads 
directly to the anomalous term in the quantum algebra. One can also repro- 
duce this anomaly by introducing a regulator (either a cut-off dimensional 
regularizationjT]), and verifying that the algebra does not close even after 
the regulator is removed. 

In the present paper, we will use a regularization that seems appropriate 
in the quantum gravity context, namely Bohr quantization 

5 Bohr Quantization 

We now quantize the partially gauge-fixed theory in the discrete Bohr quan- 
tization framework jSJ |2]. We first smear the observable X{x) with a test 
function f{x): 



Secondly we define a basis |ai, 02, a^v > of orthonormal states, such that: 



on which we define the action of the quantum smeared operator Xj: 



where L is a constant of dimension length which will be determined below. 
The discretization {xi} is at this stage completely arbitrary. The are the 
eigenvalues of X at the point Xi. This spectrum is, by construction, discrete, 
and will be determined once we choose a representation for its conjugate 
momentum, P. 

Instead of quantizing P directly we first of all note that it is the di- 
mensionless variable GIP that appears in the Hamiltonian constraint. We 
therefore define: 




(27) 





Xf\ai,a2, ....,aN >= L^(/(xi)ai)|ai, 02, > . 



(29) 




(30) 



where the constant 7 is dimensionless. The Poisson bracket is: 




(31) 
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We define the quantum version of U^{x) by its action on tlie basis states: 

U^{xi)\ai,a2, ....,a7v >= |ai, ...,ai - 7, ....,aN > ■ (32) 
Tlie definitions (j29|) and (j32|) imply the quantum commutator: 

[Xf,U^ix,)] = -^LU^ix,). (33) 

A comparison of (j31|) and (j33|) determines the arbitrary constant L to be 
proportional to the Planck scale: 

L = hGl. (34) 

We are now in a position to determine the spectrum of X. In particular, 
given the action of JJ^{xi), one has that: 

|ai, a,, Cat >= {(j±.^{xi))'^\ai^ ...^af \ ...,aN >= \ai, ...,af'^±m'y, ...,a]\i > ■ 

(35) 

Given the physical interpretation of X in terms of the square root of the 
conformal mode of the metric, it is reasonable to assume that the spectrum 
is symmetric about 0, and take a-°^ = 0. Henceforth we assume that = m7 
where m is any integer. 

We implement the Hamiltonian constraint using the following smeared 
operator: 

dfX^\au....,aN> = ^Jdxf{x)[fG^P^X^-j{<j)f-e{<j))X^)\au....,aN> 
= ^^fixi)(a^G'^l'^P^~ji-eia'^)\ai,....,aN>- (36) 

where for ease of notation we henceforth denote ji := j{(l){xi)) and similarly 
for ei. 

We now need to define P'^{xi). Following Ashtekar et. al. ^: 
/^P^(xi)|ai, ...,aAr > = Yim.^{2-U^{xi) -U^^{xi))\ai,...,aN > 

= lini ^ (2|ai, > -|ai, ctj - 7, «iv > 
7^0 7^ 

- |ai,...,ai + 7,aAr >) . (37) 
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We are now ready to solve the Hamiltonian constraint 

dfX^\ip >= 0, (38) 
where \ip > an arbitrary state in the Hilbert space: 

lip >= 'ipiai, ■■■aN)\ai, ...,aN > . (39) 



dl ,. . . ^Ctn 



Since there is no couphng between spatial points in the Hamiltonian con- 
straint, it is reasonable to separate variables: 

ip{ai, ..mn) = V^i(ai)^2(a2)---V'7v(a7v). (40) 

By substituting and (jlOI) into (jH^ and using the definitions pHjl and 
(jnZj) one finds that the coefficients ipi{ai) must each satisfy: 

/ 2 \ 

-I {^i{ai + 7) + i^i{ai - 7)) = ( - ^i^^i ~ fi\ ^i(a»)- (41) 

A normalizable state > satisfies 

<^>2:=^|^-|2 <oo. (42) 

m 

We henceforth suppress the index i labelling the spatial point. As in |3] we 

use the fact that = m'-f to write ip{ai) = tp"^ so that ip{ai ± 7) = ip'^^^. 
Hence (14 ip can be written as: 

(m7)2 [^"^+1 + ^'"-i - (2 - £7^) + 72J = 0. (43) 

To best understand the issue of convergence, we follow the lead of Ashtekar 
et. al. and Fourier transform in m to a continuous variable k, so that the 
difference equation becomes a second order differential equation. 

Multiplying ()43|) by exp{ikm'y) and summing over m, and using 

^(A;) := ^^™e*'"'=^, (44) 

m 

so that 

5^(m7)2^"^-ie'™'^^ = -dl (e-"'^^P{k)) , (45) 
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and similarly for ip"^^^, we get 
dl 



7'\ 

1 — cos (^7) — e— 1 ip{k) 



2 '2 

+ = 0. (46) 



Now choose a new variable z := 'yk/2. Then the above can be written as 
(1 + A cos 2z)i;"{z) + B sin 2zij'{z) + {C + D cos 2z)i;{z) = 0, (47) 
where the 'constants'- that is, parameters independent of z, are defined by 



A 

B = D 
C 



-(l-e(0)7V2)-\ (48) 
4(1 - e(0)7V2)-\ (49) 
2(j(0))'(l-6(0)7V2)-^ (50) 



The differential equation Eq. ()47|) is recognized as Ince's differential equation 

It is proved in ^18^ that |y4| < 1 is a sufficient condition that there are an 
infinite number of solutions of Ince's equation which are periodic with period 
71 or 27r. By Parseval's Theorem, it follows that these periodic solutions of 
Eq. ()47|) generate normalizable states satisfying the Hamiltonian constraint 
Eq. (jHS|) . The condition \A\ < 1 is satisfied for j > 2GMI, i.e. exterior to the 
horizon and when 

j < 2GMI - 4/7^ (51) 

This implies that we are guaranteed the existence of normalizable solutions 
everywhere exterior to the horizon. However, the width of the region inside 
the horizon where solutions are not guaranteed is, in terms of the area (pu = 
j-^{2GMl) of the horizon,: 

^ (52) 



(t)H 7^ 2GMI ■ 
This is narrow for macroscopic black holes with mass: 

4 

2GMI > — , (53) 

T 

^In [21 the non-relativistic simple harmonic oscillator is quantized in the above manner, 
and the differential equation corresponding to Eq. (|46(l is the Mathicu equation. 
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but is of order one for microscopic black holes. In the classical limit, 7 — 0, 
the entire interior of the black hole is potentially excluded. 

In the next section we will discuss in some detail the spectrum exterior to 
the horizon of the solutions of the difference equation. First we will present a 
heuristic argument that one expects a classical spacetime geometry to emerge 
in the region far from the event horizon. 

Consider the factor a := 1 — e(0)7^/2 which occurs in the coefficients of 
Eg . fl47p . We wish to consider the region for which a >> 1, which in turn 
implies that: 

(j-2GM/)^»l, (54) 
In this limit the solutions should approximately obey the differential equation 

i)" + - sin i2zW + AV = 0, (55) 
a 



where A := \j2/ai{(j)) which is large for large j(0), so that one can drop the 
term in -0'. In this case, the solutions are of the form 

iIj[z) = ipQsm^Xz). (56) 

and the spectrum is 

A = 2iV, (57) 

for some integer A^. Using the definition of a this gives the following equation 
for j in terms of A^: 

/ = 2iV2(^l-(2GM/-j)y). (58) 

Under the assumptions made above, the term proportional to j dominates 
the right hand side, so that the approximate spectrum far from the black 
hole is: 

3 = NW- (59) 
Using dSni) and the fact that j » 2GMI, we find: 

- = If « 2^^, (60) 

which is small for macroscopic black holes as expected. 
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6 Details of the Spectrum 



We now establish an asymptotic expansion for each of the infinitely many 
allowed values of the spatial coordinate as parametrized by the quantity 
We consider the difference equations Eq. pH|) (or. equivalently, the re- 
cursion relations for periodic solutions of Ince's equation). These equations 
are rewritten as 

= 0, (61) 
^n+i^^n-i ^ (2a--^)7/>", n=l,2,... (62) 

posed as a second order difference equation and one initial condition. In the 
above, 

a := l + (j(0)-2G'M£)7V2; 

b := j2(0). (63) 

We will obtain an infinite number of self-consistency conditions, fixing the 
spectrum. 

Since the above is homogeneous and linear in the tp^, we may normalize 
by taking tp^ = 1 and then obtain a general solution by multiplying Eq. ()62p 
by an arbitrary constant. With these two initial conditions the tp^ are given 
for all nonnegative integers n. We require further that the ip" converge for 
large n in order to assert that an infinite set of equations is satisfied. The 
convergence condition is the constraint on a and b that we are seeking. When 
we have this condition the ip"' will be dictated by conditions near n = oo and 
the difference equation then gives all the ip"^. When ip"^ becomes known, then 
use of the initial conditions and the difference equation for n = 1 gives the 
constraint ip'^ = 2a — b. The order of the difference equation is lowered by 
noting that it involves only the expression ip"'~^^/ip"' = kn and making this 
substitution leads to the first order difference equation and initial condition 

kn = 2a — b/n^ — - , n = 2,..., 

kn-l 

ki = 2a -b. (64) 
The constraint on a and b is now given by /(a, b) = where 

f{a,b) = 2a-b- ki{a,b). (65) 
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From Eq. (j64|) we have 



2a — h/n"^ — kn 



so that the constraint ()(i5|l can be written as a continued fraction equation: 

= V, - (67) 

where: 

Vn:=2a--^. (68) 

Let us now consider the case that a and b are large, both because this 
condition leads quickly to the algorithm for finding the relation between a 
and b and because the numerical expression of this relation suggests what 
asymptotic expansion to look for. If n is very large then the term b/n"^ can be 
neglected and we find that Eq. ()66|) is, in this limit of large n, a very quickly 
convergent iteration. If we are to solve k = g{k) by iteration, then a fixed 
point is convergent if g'{k) is less than one in size and for the above iteration 
we have rapid convergence. The algorithm for evaluation of /(a, b) is to find 
kn-i for large n, given a and b and repeatedly substitute until ki is found. 
Since 2a — b/n^ is large and kn is small we start the iteration with the guess 
that kn = 0. At the same time, since one generates the sequence of values of 
kn, one can obtain a sequence of derivatives 

, kn,a 2 kn^b + l/n? 

i^n—l,a T~f o r To? '^n—l,i 



{2a - b/n^ - A;„)2 ' (2a - b/n'^ - knf ' 

again starting the iteration with the derivatives at n being zero. Newton's 
method can then be performed on the equation /(a, 6) = and the result 
is plotted in Figure 1. This diagram, it may be noted, has curves that for 
large a and b appear to straighten out and their slopes happen to be 2m? for 
integers m = 1, 2, . . .. 

We now seek an asymptotic expansion for the curves of Figure 1 for large 
a and b. We notice that values for kn propagate stably from high to low 
values of n by substitution into Eq. (|M|l since all kn are small and the value 
of kn-i is dominated by the a and b terms of the denominator. On the other 
hand, we note from Eq. ()64|) that unless the large a and b terms happen to 
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cause the difference 2a — h/ri^ to be small for some integer n, first, ki is large 
due to the initial condition and then /c„ is large for increasing integers n. 
Thus, we obtain a contradiction unless h is close to 2am? for some integer m. 
The algorithm for solving for the fc„ is to find the small kn for n > m and 
the large /c„ for n = 1, 2, ... m. The coefficients for km found both ways must 
be consistent and their equality gives the conditions that determine the 6j. 

The following substitution happens to give the correct form for the solu- 
tion to the difference equation in Eq.(|Mj). Set 

oo oo 

fcn = Ec.(n)«'"'\ b-2an^ = Y.ka'-^\ (69) 

i=0 i=0 

where bo = 2{m? — n^) and the rest of the hi are independent of n and 
unknown. First we determine the form of the A;„ that propagate in from 
large values of n. Since fc„ is small we know beforehand that Co(n) = for 
n>m. Placing the form Eq.(|69|) into Eq. (j66j) . we have 

oo r 

{ci{n)cr-i{n - 1) + Cr-i{n - l)hi/n^^ a^"^'" + 6i^r = 0. 

r=0 i=0 

Equating coefficients of each power of a we obtain 

ci(n - 1) = — ^ , n = m + 1, m + 2, . . . . 

2(1 — m^/ra"^) 

Ci{n)cr-iin-l) + Cr-i{n-l)bi/n'^ ,i n 

Cr{n-1) = 2(1 -mVn^) ' ^ = 2'3'---' n = m+l,m+2, 

(70) 

We now propagate solutions for kn beginning at n = 1. Substituting Eq.(|U^ 
into Eq.(|Mj) we obtain the result 

^ ^ _ ELo Cr-i(n^l)b,/n^ - 6r,l + E[=0 Ci{n)Cr-i{n - 1) _ . ^ 

Cr[n)- co(n-l) ' ^-U'^'^'--- 

(71) 

which may be used for n = 2,3, ... ,m to give a second expression for the 
coefficients of fc„ in addition to that of Eq. ()7()|l . Equating the results of 
Eq. (f7n|) and Eq. (f7T|) allows one to solve for the parameters 6j. The leading 
examples follow. 

Table 1 The Asymptotic Expansion of 6 as a function of a. 
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m relation between a and b 



1 h-2n - ^ i2 889_ 

1 u — z,u 3^ ^gg^3 9720a5 



2 6 = 8a - 

3 6 = 18a 



44 2626 2157091 



15a 3375a3 5315625a5 
657 27080019 



4 6 = 32a- ^, 

baa 

5 5 = 50a - 

231a 



70a 10976000a» 
1310 



1751309 



6 b = 72a — 

27456a 

The above establishes Eq. ()58|) to lowest order, and gives the first order cor- 
rections. Moreover, the definitions of a and b imply the relationship: 



( 



2 Y 

2GM/ + — (a-l)j . (72) 



For fixed M and 7,the expansions in Table 1 in conjunction with the 
relationship (|72|l can be used to determine the allowed spectrum for ji, and 
consequent ally 0(xj), i.e. the discretization of the spatial slice. For large 
Table 1 confirms the approximate spectrum (j^HI), which approaches the 
continuum for large quantum numbers m. For small and hence a close 
to 1, one needs to keep a suitable number of terms in the corresponding 
asymptotic expansions. 

Figure 1 shows qualitatively the effects of changing the mass M of the 
black hole and the discretization scale 7. Increasing the mass shifts the 
parabola down (small to large dots) and the spectrum for a is squeezed 
closer to the horizon (a = 1). Decreasing 7 (dot-dash to large dots) flattens 
the parabola and moves the eigenvalues of 6 = (j(0))^ closer together. This 
behaviour corresponds to our physical expectation that the spatial topology 
near the horizon becomes more continuous with increasing mass or decreasing 

7- 

Note that the mass is expressed in terms of Planck units. For very large 
M, one expects the spectrum in the exterior of the black hole to be approxi- 
mately continuous even near the horizon. This in fact can be seen to happen 
by noting from (j72j) that as MGl increases, the horizon a = 1 occurs for 
larger and larger b, which in turn implies that the parabola defined by (f7^ 
first intersects the asymptotic expressions in Figure 1 for larger and larger 
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m, where the hnes are more and more dense. Thus the lowest eigenvalues 
get closer and closer together as MGl increases. 

7 Conclusion: The Quantum State of a Black 
Hole 

We have presented a quantization scheme for the spacetime surrounding a 
generic, single horizon, spherically symmetric black hole. The key features 
of the scheme are: 

• Only the spatial diffeomorphisms were fixed at the classical level; the 
choice of gauge fixing allows classical slices that are non-singular across 
the future horizon. 

• The single remaining constraint (the Hamiltonian constraint) docs not 
contain any spatial derivates. It describes the dynamics of a particle 
moving in a ^ potential, where the role of particle position is played 
by X{x) — ef, the square root of the conformal mode of the metric. 

• The Hamiltonian constraint is quantized using Bohr quantization, which 
by construction yields a discrete spectrum for X and induces a discrete 
spatial topology. 

The interesting result is the induced discretization of the spatial topol- 
ogy exterior to the horizon. It is important to remember that our procedure 
specified the spectrum of the operator X at each spatial point. The Hamilto- 
nian constraint then determined the spectrum for jj, i.e. the discrete spatial 
topology. Moreover, this topology has the desirable property that it approx- 
imates the continuum far from the horizon even for microscopic black holes, 
as well as near the horizon for macroscopic black holes. 

We end on a speculative note. Prom Figure 1 it is clear that the parabo- 
las intersect each line twice; the intersection furthest from a = 1 determined 
the spectrum of to which we referred to above. However, the second in- 
tersection, closer to the horizon may also have physical significance. Note 
in particular that each parabola (i.e. for each value of M and 7, there are 
an infinite number of intersections infinitesmally close to the horizon. This 
seems to imply that although the spatial topology just outside the horizon 
is discrete, there is in fact a continuum field theory residing on the horizon 
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itself. At each spatial point in this (approximate) continuum theory, there 
is a set of SO (2,1) generators that are defined classically in terms of the 
operators X and P. It may be that the anomaly associated with the quan- 
tum version of this algebra near the horizon is connected to the statistical 
mechanical entropy of the black hole, much as in the work of [0] and |lUj . 

There is of course a great deal yet to be done. The spectrum of X and 
the resulting spatial topology below the horizon must be understood and the 
properties of the wave function at each spatial point should be analyzed. A 
complementary approach to the quantization of the same system, namely 
Schrodinger quantization is also being explored [T^ and this should also yield 
insights into the physical interpretation. Most importantly, it is of great 
interest to see whether the present formalism can be implemented in the 
presence of matter, which would provide a new mechanism for analyzing the 
quantum dynamics of black hole formation and Hawking radiation. All this 
and more are currently under invest igation|2Uj. 
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b 

Figure 1: Solid lines graph the relations between a and b from asymptotic 
expansions in Table 1. The parabolas graph (f?^ for various M and 1/7^. 
The large dotted line is 2GMI = 1.90, 1/7^ = 1.77; the small dotted line is 
2GMI = 1.00, 1/72 = 1.77; and finally the dot-dash line is 2GMI = 1.90, 
1/7^ = 1.18. The horizon is at a = 1. Intersection points of parabolas with 
the solid lines yield the spectrum for either in terms of a := l + ^(j(0) — 

2GM£) or of 6= 
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